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ABSTRACT 


A model has been developed for estimating the number, 
sizes, and interfluvial distances of streams to be crossed 
in traversing a naturally occurring drainage basin. A 
relationship among four network parameters has been found 
which, upon quantification, provide sufficient information 
to completely structure a mature drainage system. The four 
parameters are: stream length ratio, bifurcation ratio, 
basin shape parameter, and coefficient of drainage density. 
The model is amenable to computer solution. Procedures are 
described for parameter quantification using common topo- 


graphic maps. 
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lee SUMMARY 


In order to determine requirements for military stream 
crossing equipment and materiel, it is necessary to deter- 
Mine the number and sizes of streams (or gaps) that have 
to be crossed. The task of actually counting and measuring 
streams throughout large areas of the world either by 
reconnaissance, from topographic maps, or by other means 
such as aerial photographs, is fraught with difficulties of 
time, money, and access. Consequently, a model is needed 
that will circumvent these difficulties and provide reliable 
data on stream size and frequency in a manner sufficiently 
timely to facilitate tactical and/or strategic military 
planning. 

This thesis presents such a model, based upon previous 
theoretical and empirical work in the field of quantitative 
geomorphology. A number of relationships are developed 
that were not heretofore available. The model combines two 
methods of stream ordering to describe a drainage basin in 
terms of its component streams and basins. The model relies 
on quantification of four parameters to estimate the fre- 
quency and sizes of streams within a naturally occurring 
drainage area. The four parameters are: stream ratio (a), 
bifurcation ratio (8), basin shape parameter (h), and 
coefficient of drainage density (k). 

Reliability of the model is materially enhanced by 


quantifying all four parameters within the geographical 


area of intended model use. However, the model reliability 
is little reduced if h and k are assumed to have equilibriun 
values and only a and 8 are quantified within the area of 
intended model use. Consequently, employment of the model 
would be greatly facilitated if topographic maps of the 
appropriate areas were overprinted with the applicable a 

and 8 values. 

The model does not attempt to predict the hydrologic 
characteristics of streams such as depth, velocity, width, 
etc., for the different stream orders. Such a model is a 
desirable adjunct to the model presented by this paper and 
offers fertile ground for further research. A model of 
this general nature that could be used was formulated by 
Mayer [Ref. 32], however, further study needs to be done 
regarding the essential features of bank conditions and 
the effects of the works of man such as bridges, lock and 
dam restrictions, etc. It is believed, nonetheless, nee 
the model presented by this paper is an essential step in 
obviating detailed ground, map, and aerial reconnaissance 


for large-scale planning and capabilities studies. 


Li. INDRODECTION 


With the possible exception of the enemy, perhaps the 
greatest single obstacle to cross-country movement of a 
combat force has traditionally been the naturally occurring 
stream. In fact, the location and obstacle characteristics 
of streams have on more than one occasion practically dic- 
tated the nature of the tactical plan. But, with the state 
of modern stream crossing technology it is more instructive 
to think of such obstacles in terms of the relative degree 
to which they impede the speed of movement of the combat 
force. In this context the problem assumes more general 
characteristics and has been properly labelled the "gap 
crossing problem." 

A gap is defined as any cavity in the earth's surface 
caused by flowing water. Thus, the study of gaps is neces- 
sarily the study of streams. Streams are of many shapes 
and forms but may be generally classed as ephemeral, inter- 
mittent, or perennial. An ephemeral stream flows only 
during or immediately following periods of precipitation 
and handles only surface runoff. An intermittent stream 
flows only during the wet season and differs little from 
ephemeral streams. A perennial stream, on the other hand, 
flows constantly. The perennial stream is generally con- 
sidered to have the greatest obstacle value. However, with 
modern stream crossing means, other aspects of the streams 


are often as important as the presence or absence of water. 


In any event, streams (or gaps) have sufficient ob- 
stacle value that no tactical plan involving the cross- 
country movement of troops is complete without a detailed 
consideration of the streams within the zone of movement. 
This has traditionally been accomplished by examining the 
area to be traversed by personal reconnaissance, examining 
topographic maps, or, more recently, aerial photographs or 
electronic means. However, it is not difficult to visualize 
circumstances under which neither method will produce 
satisfactory results. For example, personal reconnaissance 
in prospective enemy territory during peacetime is quite 
difficult in addition to being quite limited in scope. 
Further, topographic maps of the tactical variety seldom, 
if ever, provide sufficient detail for the smaller gaps. 
Aerial photographs and the other, electronic, means 
suffer from these shortcomings plus being gquite expensive 
and time-consuming to obtain and interpret. In short, a 
better method is needed for determining (or estimating) 
the number and sizes of gaps within a given area, a method 
that is inexpensive, timely, and widely applicable. This 
model, then, seeks to provide that method by describing the 
behavior of naturally occurring stream networks in mathe- 
matical terms and using the theory of the most probable 
state to secure useful predictive results. 

The study of the formation of drainage patterns from a 
stochastic viewpoint had an early beginning. In 1802, John 


Playfair, Professor of Natural Philosophy, University of 


Edimiburgh, wrote Tllustrations of the Huttoenian Theory of 
the Earth. The following passage has been called Playfair's 


Law: 
Every river appears to consist of a main trunk, fed 
from a variety of branches, each running in a valley 
proportioned to its size, and all of them together 
forming a system of valleys, communicating with one 
another, and having such a nice adjustment of their 
declivities that none of them join the principal valley 
either on too high or too low a level; a circumstance 
which would be infinitely improbable if each of these 
valleys were not the work of the stream which flows 
cbt. 
Although Playfair obviously recognized the systematic aspects 
of a natural drainage system, little was done to quantify 
these aspects prior to the work of R. E. Horton in 1945 
[Ref. 11]. Horton's bifurcation and stream length ratios 
are particularly important to this model and to most of 
the work that has been accomplished in the field since 
Horton's publication. In recent years probability theory 
has been used in an attempt to explain and quantify various 
aspects of naturally occurring drainage networks. Prominent 
among the works in this area have been those by Hack [Ref. 
27], Leopold and Langbein [Ref. 28], and Schenck [Ref. 18]. 
The general theory basic to these works is that the forma- 
tion of drainage networks can be explained by the laws of 
probability and that, free of constraints, streams occur in 
a random fashion. Leopold and Langbein viewed the laws 


of Horton as representing the most probable state in a 


stochastically formed network, and used a random walk 


model to construct a network representing the mest probable 
state to show that the laws of Horton did, in fact, hold. 
Schenck used a similar model to substantiate the empirical 
work of Gray [Ref. 9]. The first attempt at constructing 

a predictive model was made by Hugo E. Mayer [Ref. 32]. The 
Mayer model was an attempt to use the work of the several 
previous authors in the field to construct drainage basin 
relationships that, with quantification of a few parameters, 
could be used to predict the numbers, orders, and hydrologic 
characteristics of streams within a naturally occurring 
drainage network. The model presented by this thesis has 
some of the characteristics of the Mayer model, but is 
designed to be a more sophisticated and a more exact 
presentation. 

This thesis will seek to build upon accepted laws and, 
using these laws and other considerations, to establish a 
number of relationships not heretofore available. The laws 
of Horton [Ref. 11] and Schumm [Ref. 34] will be used 
essentially as presented by these authors with an additional 
law to be deduced from the empirical work of Hack [Ref. 27]. 
The first effort toward model construction will be to define 
an admissable region for the key model parameters. The 
polynomials describing the upper and lower bounds for the 
region were defined by Mayer [Ref. 32]. The polynomial 
describing equilibrium behavior within the region will be 


developed by making certain assumptions concerning the 
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"Jumping" process and then describing in mathematical terms 
the proportion of streams of a given order that jump and, 
of those that jump, the proportion that jumps to a specific 
order. The product of the latter two proportions produces 
the proportion of streams of order i that jump to order j 
fow all i = 1, 2, ---, N, and all j > i. In particular, 

it will facilitate describing the proportion of streams 
within a drainage area of order i that jumps to order N, 
the order of the parent stream. This relationship employed 
in conjunction with an equation by Hack [Ref. 27] will be 
used to write the polynomials describing the admissable 
region for the model parameters. This admissable region 
will be used to estimate parameter values within drainage 
@eeas for which specific parameter quantification is not 
practical and to estimate the error to be expected from such 
estimates. In defining the admissible region a number of 
assumptions will be necessary. The thesis will attempt to 
evaluate the impact of these assumptions on ultimate model 
output and to establish relatively simple relationships 
that demonstrate the logic underlying the assumptions. 
Finally, the model will be given a statistical structure by 
assuming distribution functions for the relevant variables. 
A uniform stream distribution within the drainage basin will 
be used to produce a sample solution, however, with rela- 
tively minor modifications the model will accommodate other 
distribution functions. Mean conditions will be used to 


write equations describing the various aspects of importance 


ila 


to this model and the confidence to be associated with each 
estimate. A five-percent perturbation will be used to 
demonstrate the sensitivity of the model to its various 
parameters and to establish those parameters for which 
Specific area quantification is necessary for a high measure 
of model reliability. The latter will facilitate recom- 
mendations concerning exogenous preparations required for 


effective model employment. 
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TII. BASIC RELATIONSHIPS 


NOTATION 

= stream length ratio 

= bifurcation ratio 

= overland flow ratio 

= basin area ratio 

= basin shape parameter 

= drainage density 

= coefficient of drainage density 

= order of the main stream of the parent basin 
= an index denoting stream order (i = 1,2,---,N) 
= an index where j} = N - i (3 = 0,1,---,N-1) 

= proportion of ith - order streams that jump 


= preopomtieon“of ith - order streams that jump to 
SeLoams “eLrmonrder 


= mean area of ith - order basins, square miles 

= area of the kth basin of ith order, square miles 

= total area of all ith - order basins, square miles 
= number of ith - order streams 

= length of the kth stream of ith order, miles 

= mean length of all ith - order streams, miles 

= length of the parent basin, miles 


= mean interfluvial distance, i.e., distance between 
streams of ith order, miles 


= width of the parent basin, miles 


= mean area drained by overland flow occurring directly 
into an ith - order stream, square miles 


is 


FEF. = total area drained by overland flow occ@ummring direc, 


into an ith - order stream, square miles 

n¥ = the number of ith - order streams encountered by a 
force in one traverse of a mean basin along the minor 
axis 

k 

i the number of ith - order streams encountered by a 
Force in one traverse of a mean basin along the major 
axis 

* 

o. = expected total length of ith - order streams encoun- 
tered by a force in one traverse of a mean basin 
along the minor axis, miles 

kx 
d. = expected total length of ith - order streams encoun- 


tered by a force in one traverse of a mean basin 
along the major axis, miles 
B. DEFINFIONS 
1. Drainage Basin. The area of slopes which lead into 
a Single stream system 
2. Drainage Divide. The line of separation that 
divides the precipitation that falls on two adjoining basins 
and directs the ensuing runoff into one or the other river 
system. The definition of drainage divide used in this 
thesis will refer to the topographic or surface divide as 
distinguished from the phreatic or underground divide 
3. Stream Order. The stream order is a measure of 
the stream's rank in the hierarchy of streams that consti- 
tute the drainage network 
4. Basin Order. The hierarchal order of a basin is 
the same order as the highest order stream in the basin 
5. Parent Basin. The parent basin is the drainage 


basin containing the highest order stream, and all lower 
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order tributaries with their associated basins, of the 


drainage system being studied 


C. STREAM ORDERING METHODS 

In developing this model two methods of stream ordering 
will be used. The one presented by Strahler [Ref. 23] will 
be used to develop the bifurcation ratio and the one pre- 
sented by Horton [Ref. 11] will be used to develop the 
stream length ratio. Both methods are illustrated in 


Pagure 1-Tf. 


HORTON STRAHLER 


3 —_— 





Figure 1-1. Horton and Strahler definitions of stream order 
applied to channel network of Hightower Creek, upper 
HiwasSee River, Towns County, northern Georgia (ex- 
tracted from Shreve [Ref. 19]). Order indicated by 
number near upstream end of respective streams. Unnum- 
bered streams are first order. If streams A, B, C, or 
D were actually second order, rather than first, then 
network would be fifth order, rather than fourth. 


ins 


Comparison of the two methods as illustrated facilitates a 
rapid aie forenceimen of the results provided by the 
different stream ordering methods. The Horton method traces 
each stream to its drainage divide and considers that each 
stream extends from its mouth to primary headwater. Thus, 
it is intuitively appealing to associate the stream length 
ratio with this method. The Strahler method, on the other 
hand, defines a first order stream as one having no tribu- 
tarles and provides that when two first order streams 
intersect, a second order stream is formed, that the inter- 
section of two second order streams forms a third order 
stream, and so on. Thus, the streams are not traced back 
to primary headwater but only to the intersection of two 
streams of the next lower order. The latter system proves 


quite useful in developing the bifurcation ratio. 


De BASIC LAWS OF DRAINAGE SYSTEM EVOLUTION 

The model will draw heavily on two basic laws presented 
by ‘Horton (Ref. 11], one by Schumm [Ref. 34], and another 
derived from works by Hack [Ref. 27]. The Horton laws are: 


(1) Law of Stream Numbers: The numbers of streams of 
different orders in a given drainage basin tend closely 
to approximate an inverse geometric series in which the 
first term is unity and the ratio is the bifurcation 
Bacto 


(2) Law of Stream Lengths: The average length of 
streams of each of the different orders in a drainage 
basin tend closely to approximate a direct geometric 
series in which the first term is the average length 
of streams of the first order [Ref. 11:291]. 


The Schumm law of stream areas states: 
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the mean drainage basin areas of streams of each order 

tend closely to approximate a direct geometric series 

in which the first term is the mean area of the first 

order basins [Ref. 34:14]. 
Hack observed that overland flow contributes some runoff 
directly into streams of all orders. The length of runoff 
directly into streams of order greater than one is, of 
course, insufficient to sustain a first order stream. 
Otherwise, a first order stream would be formed. However, 
Hack observed that the area drained in this manner was 
proportional to the stream length ratio. Hack's observa- 
tion implies our final law of drainage evolution which can 
be expressed as: 

the average areas drained by overland flow directly 

into streams of each of the different orders within a 

parent basin tend closely to approximate a direct 

geometric sequence in which the first term is the mean 

area of the first order basins. 
BE. THE MODEL RATIOS 

The foregoing laws enable us to define the following 
dimensionless relationships. 

1. Stream Length Ratio. The stream length ratio is 


the ratio of the average length of streams of a given order 


to the average length of streams of the next lower order: 





ye, 
a= — ee (11 
on 
2 Eitumeeitien Ratre. The bifurcation ratio is the 


ratio of the number of streams of a given order to the num- 


ber of streams of the next higher order: 


lee 


— 
B = ee a (1 


Nis] 





3. Overland Flow Ratio. The overland flow ratio is the 
ratio of total area drained by overland flow directly into 
the streams of a given order to the total area drained by 
overland flow directly into the streams of the next lower 


OEGen: 





a (1.3) 


4. Basin Area Ratio. The basin area ratio is the ratio 
of the total area of basins of a given order to the total 


area of the basins of the next higher order: 


i= 2 ae (1.4) 





A hypothetical third order basin using the Strahler method 


of stream ordering is illustrated in Figure 1-2. 


Fes A DESCRIPTION OF THE. BAST Cerri 


1. Horton's Laws. Horton's’ £irst law.can be wrueten as 
. Nel 
a 2 N-1 
{Bo} = By > - : «cee ’ 
i=0 


where the first term is the number of streams of order N and 
the last term is the number of streams of order l. Thus, 


the number of streams of order P ean be wWwrreeen acs 
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EL fl 2 


First order basin not included in a parent 
second order basin 


First order basin included in a parent second 
order basin 


Second order basin portion not included in a first 
order basin (overland flow directly into second 
order stream) 


Third order basin portion not included in a first 
or second order basin (overland flow directly 
into third order stream) 


Figure 1-2. Basins of Hypothetical Third Order 


Stream (Redrawn from Mayer [Ref. 32]) 


1G 


nN. =e ! : (Les) 





we see that 8 is independent of basin order and thus can be 
considered constant for all i< N. Horton's second law can 


be written as 


where the first term is the average length of all first 
order streams and the last term the length of the Nth - order 


Stream. Thus, the mean length of an ith - order stream is 


given by 
Q. = Zar"? (1.6) 
where 
ny 
a. = wily > Lay 
Ni k=l 


If we form the stream length ratio 
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we see that a is also independent of basin order and can 
be consicaered constant for all 1 > 1. 
2. Schumm's Law. Schumm's law of stream areas can be 


written symbolically as 


. N-l 
(2,0 diye = 2p aoe el ,, 2. ae 
= 
wieke the first term is the mean area of all first order 


basins and the last term the area of the parent basin. Thus, 


the mean area of the ith - order basin is given by 


—- _ = ,i-1 
a a,o ; 


which can be written as 


at 7 
A; = —— A,8 
sai 
Note that 
i @ ge 
1 —— 
wa => 2b} @eaeee= ina 
al k=] | aga ke sake 


imeewe forum the basin area ratio 





we see that A is independent ot bhaSin order and thus con- 
Seomrereogeall i < N. By solving recursively we obtain the 


more useful relationship 


Ea 
Ay-5 area (Tae 
3. The Law of Overland Flow Areas. The law of overlane 


flow areas can be written symbolically as 


where the first term is the mean area of all first order 
streams and the last term the area drained by overland flow 
occurring directly into the Nth - order stream. Thus, the 
area drained by overland flow directly into an ith - order 


stream can be written as 


If we further define FE. as the total area drained by over- 


land flow directly into an ith - order stream, we can write 


Bo Seal ,oel : 
a ll 


Forming the overland flow ratio we obtain 


= ol 
F. n,a,o 7 al 





i 
| 
Q 
WD 
i 


Dy 


which is constant within a relatively homogeneous drainage 
basin since by Horton's first and second laws, a and 8 are 


constants. Solving recursively we obtain 


Ze 


and since F, = A., we further obtain 


Finally, wsing equation 1.7 with i = N-j, we obtain the use- 


ful relationship 


ae =i A, (1.8) 

G. CONVERSION BETWEEN STREAM ORDERING METHODS 

Since the Horton method of stream ordering extends the 
streams back to primary headwater, we must convert the 
stream lengths to the Strahler method in order to make them 
compatible with the bifurcation ratio. An examination of 
Figure 1-1 reveals that to accomplish this conversion 
requires only that we subtract from each Horton mean stream 
length of a given order a Horton mean stream length of the 
next lower order. Symbolically, if we define he = T,a0*7* 


using the Horton method, then the corresponding mean stream 


length using the Strahler method can be written as 


ee 2, (a = 6) ieee]. (1.9) 
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which reveals that the stream length ratio is unchanged by 
the conversion and, thus, the two stream ordering methods 
can be used interchangeably as far as the stream length 
ratio is concerned. The numbers of streams of various 
Orders can also be quickly converted from one method to 


the other. If we denote by n. 


4H the number of streams of 


order j uSing the Horton method and by n.. the number of 


1S 


streams of order j using the Strahler system, then 


and, conversely, 


nN. = 


5H Nas ~ N541,8 Pees ee ep Nie 


Thus, as indicated by Shreve [Ref. 19], the number of first 
order Strahler streams iS equal to the total number of 
Horton streams within the parent basin. However, the bifur- 
cation ratio is unchanged and, thus, the two methods are 
again interchangeable with respect to the corresponding 


dimensionless ratio. 
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IV. FUNDAMENTAL PROCESSES 


A. JUMPING 

We hypothesize that within an Nth - order parent basin 
streams of order i bifurcate with streams of order itl, 
1+2, ---, N, in direct proportion to the relative total 
stream length of the recipient streams. Note that a stream 
of order i cannot bifurcate with any order less than itl; 
therefore, relative stream length refers to streams of 
orders i+tl and greater. This hypothesis is supported by 
considerations of probability in a relatively homogeneous 
drainage basin. As indicated by Leopold and Langbein 
[Ref. 28] in their concept of entropy in landscape evolution, 
rainfall occurring on a uniformly sloping surface is as 
likely to produce a stream of a given order on one part of 
the surface as on another and the orders of bifurcating 
streams are largely determined by chance. Thus, it 
follows that first order streams, for example, will inter- 
sect third order streams in direct proportion to the relative 
total length of all third order streams within the drainage 


basin, where relative stream length refers only to those 


streams of order greater than one. Under the Strahler 
system, at least two ith - order streams are necessary to 
form each (i+tl)st - order stream. Therefore, we can express 


the number of ith - order streams that bifurcate with 


(i+1)st - order streams as 


Zo 


7s oe eee) 
itl N : 

b See 

j=it1 / 7 


Ne gen 7 2MG47 FT OYA 


— 


Dividing this by n, we can write the proportion of ith-order 


streams that bifurcate with (1i+l)st-order streams as 


Rey item “On - (ne-2ns 4a) Dey ea ay 
ne gle ne N j 

= 2 - Se ee. 

j=it1 J J 


Since the proportion of ith-order streams that jump is the 
complement of the proportion that bifurcates with streams 


of the next higher order, we can write the proportion that 


jumps as 
oe ee (ns-2ngig) M54. 447 
a n, n. N —— 
> ane 
j=it1 J J 


which, using equations 1.5 and 1.6, reduces to 


N-i -l 

(p-— Ns ll se) 
P. = ————-—.————— , i < N. (2, 1) 
ep 
Equation 2.1 describes the behavior of P. within certain 
bounds. The lower bound is trivial and the upper bound is 
determined from the fact that at least two ith-order streams 
are required to constitute each (it+l1)st-order stream. The 


bounds on P. can be written as 
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which can be further expressed as 


—, 
M < Ea i = 2hae a 


The behavior of P. in a tenth-order basin wherein a and Bare 
arbitrarily chosen to be 2.68 and 4, respectively, is illus- 
trated in Table 2-1. The table also gives the proportion 

of streams that don't jump as qa In evaluating the table 
contents note the important characteristic that the only 
stream permitted to exit a parent basin is the Nth-order 
stream. All other streams lie entirely within the Nth~-order 


basin. 


TABLE 2 


BIFURCATION BEHAVIOR OF INCLUDED STREAMS WITHIN 
A HYPOTHETICAL TENTH-ORDER DRAINAGE BASIN 


ORDER PREEORTION THAT PROPORTION THAT 


(i) JUMP (P.) DON'T JUMP (q, ) 
1 . 336 664 
2 Fece 667 
3 .329 SO 
4 as: 677 
5 oud . 686 
6 aa = “iG 
7 267 Be: 
8 ite a 
9 0 1.000 


a7 


Now, of the ith-order streams that jump, we define Pan as 
the proportion that jumps specifically to order r. The 
process of jumping within a parent basin is then as illus- 
trated in Figure 2-l. Further, we can write, symbolically, 


the proportions that jump to order N as 





en 
_ he 
PN-2,N — n er, 
NN 
- N'N _ — 
N-3,N = — 7 -l — 
yn ee a Se 
ty Dy iT 
and PN-3 a ra =F, 2-3 
a 
NTN Sy O36 an 44+2"n-5+2 i tO, +. scone 
1 - ae 


The proportion of all streams of order i that jump to order 


N can now be written as 


es nate desea Ses er! ee 
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Be THE BASIC PARAMETRIC RELATIONSHIP 
Schumm [Ref. 34] added to his law of stream areas "It 


could be assumed that such a relationship would exist if 
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Figure 2-1. The Jumping Process 
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there were any connection between the length of a stream 
and the size of its drainage basin." Hack [Ref. 27] later 
showed that such a relationship did exist. The Hack 


relationship can be expressed as 


i aoe (2.3) 
where Re represents the mean stream length using the Horton 
stream ordering method. Hack found that k could assume 
values™in™ titegrange slime ek. < 2.5, but that a value of 
1.4 was typical. Mayer [Ref. 32] observed that h was bounded 
by 0.5 <h < 1.0, but that the typical value was 0.6. If 
we assume that within a relatively homogeneous drainage 
basin the basin shape parameter (h) and the coefficient (k) 


for included basins are constant, we can write 





Multiplying both sides by go and rearranging we obtain 


g tyt7n - lL. oe 





which reduces to the basic parametric relationship 


WAgl 


B = Aa (235 


The constancy of h and k are discussed more fully in Chapter 


Jbg is 
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C. THE FUNDAMENTAL EQUATION 
Considering both jumping and overland flow, we can now 
account for all of the area composing a parent basin. In 


symbolic form the area is 


By = Ay + Py-oPy-aewtn-2 7 °°° * FoF aa * “NU 


which, using equations 1.7, 1.8, and 2.2, becomes 


2 -l 
= (ap ) (1-26 ~) 2 
An = Ay + i= A Ay + 
AE 6) 
N-2 N-l -l, 
# i ems Yale 4 (pa) *a, 
= 


and further reduces to the fundamental equation 
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Equation 2.5 defines the admissable region for a,f pairs for 
selected basin shape and size. At this pornt itis “insttruc- 
tive to note that if there was no jumping and no overland 
flow directly into streams other than first order (or, if 
the overland flow in adjacent basins was equal in area 
drained), then the basin area ratio would be unity. However, 
overland flow directly into streams of all orders is always 
present and is geometrically related to the stream length 
ratio. Thus, although A is a measure of both overland flow 
and the process of jumping, once the order of the parent 


basin is specified and the stream length ratio determined, 


Sih 


the overland flow is fixed. Consequently, since only Pi, 
as given by equation 2.1, can vary, the bounds on P. deter- 
Mine the bounds on A. Hence, A assumes its minimum value 
when P, is a maximum, i.e., Pa (B-2776, “all Gn Panay 
assumes its maximum value when P. = 0, all i<N. Therefore, 


equation 2.5 can be solved for three values of P.. With 


P. = 0, all i<N, we obtain 


\ (oS 1 (2.6) 


which defines the upper bound on the admissable region for 


the a,8 pairs. With P. = (8-2)/8, all i<N, we obtain 


a ee 0, ui) 


(4) 2 3 pace 
which defines the lower bound on the admissable region for 
the a, pairs. With P. defined as indicated by equation 2.1 
we obtain equation 2.5 which defines the most probable 

a,8 values admissable for each N. Now, by using equation 
2.4 to identify the interdependence of the four parameters 


we can write equations 2.5, 2.6, and 2.7, respectively, as 


V1 0 oat ey oa! _ A s 
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(2.10 ) 
Solution of equations 2.8, 2.9, and 2.10 yield the admissable 
region for the a,8 pairs for each parent basin and speci- 
fied h. The equations can be solved using such techniques as 
the Newton-Raphson method of succesSive approximation 
[Ref. 3]. Figure 2-2 illustrates the solutions for N = 8 and 
h = 0.55, 0.6, and 0.65. Table 2-2 reflects selected values 
for the solution of equation 2.8, the equation representing 
the most probable state. The solution of equation 2.8 
using h = 0.6 will be referred to hereinafter as the equili- 
brium solution as it represents the most probable state 


within any Nth-order basin. 
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‘EQUILIBRIUM SOLUTION OF FUNDAMENTAL EQUATION 
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Figure 2-2. Admissable Region for a, Pairs 
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V. A PARAMETRIC EVALUATION 


A. THE AREA EQUATION 
Observing that the total area drained by an Nth-order 
stream system consists solely of area drained by overland 


flow occurring directly into streams of order 


i(i = l, 2, ---, N), we can write the area of the parent 
basin as 
™ —_ = — N-l 
Ay = n)a, + Noa, a -.., seer Nvy@ 4 
N 
= ay > ee 
g—2 


which sums to 





A= a, Se, ie. 


and further converts to the fundamental Hack equation 


[Ref. 27] 


i = 


Ay = a,8 ° oer: (3.2) 


Using the latter equation and forming the basin area ratio 


we obtain 


BS) 





However, the basin area ratio as defined by equation 1.4 


requires that all areas A in size within the parent basin 


Ne] 
be considered as the numerator and not just one such area 
as described by the Hack equation. Thus, the basin area 


ratio becomes 


A = = SS (333) 





which defines the upper bound on iA since Hack did not en- 
vision the process of jumping in his development. This 
curious omission of Hack's plus the fact that his develop- 
ment does consider overland flow occurring directly into 
streams of order greater than one implies that ratio 3.3 
differs from unity by the fraction of overland flow occur- 


ring directly into the Nth-order stream, 1.e., 


It can easily be demonstrated that FV /Ay approaches zero 
rapidly as N increases. Thus, equation 3.3 is essentially 
equivalent to equation 2.5 with P. = 0, all i<N, the latter 


implying that 
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(Ap) 


(3.4) 


The closeness of equations 3.3 and 3.4 is illustrated by the 


results contained in Table 3-1 which were achieved using 


N 


a-values of 2, 


(Figure 2-2) for each a value. 
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TABLE 3-1 


moudtien 3.3 (2) 


P=.656 


0 


0.60 


Oe 7 


p=.398 


0 


de'O'0 


1.00 


1.00 


Equation 3.4 


3, and 4, and upper bound values on 8 


COMPARATIVE VALUES FOR A USING EQUATIONS 3.3 AND 3.4 


(X) 


p= .398 


0 
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0.90 
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o.o 
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LOO 


1.00 


Since the two equations are essentially equivalent, we may 


interchange them in equation 2.5 to obtain 


Sif 
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Thus, we have demonstrated the relationship between the Hack 
equation (3.1) and equation 2.5, the fundamental equation of 


this thesis. 


B. THE COEFFICIENT OF DRADNAGE DENSITY 

The geometric laws of Horton, Schumm, and Hack from which 
were deduced the parameters q,\ and € are based on the values 
of Ly and ay which for simplicity of development most 
authors assume to be unity. In this model, however, we are 
concerned as much with absolute values as we are with the 
constancy of their adjacent order ratios. Thus, it is 
necessary that we examine the behavior of ay and Ls within 
relatively homogeneous drainage basins. To do this we 
define D to be the drainage density within a basin area 
. Measured in linear units of stream length per square unit 
of basin area (usually in miles per square mile). The 
density of streams within an area is, of course, a measure 
of how well the area is drained. A well-drained area would 
have a relatively low D-value and, conversely, a poorly 
drained area would have a relatively high D-value. The 
equation for D is 


D= rae ; 


1 niseal 


| a 


a 
an 
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Mirren, using equation 1.5, 1.6, andes.) sums to 


Dia ac (3.6) 


Thus, we see that D is constant within any basin within which 
the aforementioned laws hold. Now, the relationship between 


ay and Ly within a baSin area is defined by eq@atiomr 2.3. 


The behavior of this equation for the equilibrium values of 
k = 1.4 and h = 0.6 is illustrated in Figure 3-l. Note that 


for a value of ay ia) @seSGuare Miles dy 1s.0.7 Males and 


Dis 2.27 miles per square mile. 


DENSITY 
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Figure 3-l. Drainage Density Versus Length and 
Area of the First Order Basins 


Thus, we see that the drainage density is constant and 


uniquely determined for any given pair of values for dy and 


a In an attempt to discover some reasonable values for 


1° 
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these constants we refer to a number of empirical investi- 
gations. Ina study of 13,376 fist order streams in the 
Appalachian Plateau, Morisawa [Ref. 33] found that the mean 
stream length wasmonoiee miles, the mean basSin area was 
0.034 square miles, and the mean drainage density was 2.44 
miles per square mile. This corresponds closely to the 
equilibrium solution illustrated by Figure 3-1. In a study 
of 146 first order streams along the iit Hiwassee River, 
New Yorke, Horton [Ref. 11] found that a drainage density of 
2.06 was prevalent. However, in a study of a fifth order 
stream system at Perth Amboy, New Jersey, Schumm [Ref. 34] 
found that the drainage density was 602 miles per square 
mile. In the same Te: Schumm recorded and compared 
drainage densities in a number of other areas throughout 
the United States ranging form 4.6 to the 602 in Perth 
Amboy. Obviously, then, drainage density assumes a wide 
range of values between different areas although it may 

be constant within Bisse basin. In fact, a number of 
authors have related drainage density to a texture ratio, 
the latter defined to be a measure of the closeness or 
proximity of one channel to another. In this respect, a 
book by Leopold, Wolman, and Miller [Ref. 1] contains the 
relationship redrawn as Figure 3-2. Thus, we conclude that 
the area (Ay) and main stream length (2.0) of a parent basin 
of any given order is directly related to the drainage 
density and as a result varies substantially from one 


locale to another. Herfce, any attempt to develop basin 
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‘Figure 3-2. Relation of Drainage Density to 
Texture Ratio, which is Obtained by 
Dividing the Number of Contour Crenula- 
tions on the most Crenulated Contour by 
the Perimeter of the Basin. [Extracted 
from Leopold, Wolman, and Miller [Ref.1]] 


areaS Or main stream lengths using as a base ay = 1 square 
Minit and Ly = 1 unit is likely to produce results substan- 
tially in error. Thus, the mathematical structure of our 
model will have to be a function of the actual basin area. 
However, as indicated by equation 3.6, the drainage density 
appears to remain constant within any relatively homogeneous 
drainage area. Now, in Chapter II we assumed in the develop- 
ment of equation 2.4 that the coefficient (k) in equation 
2.3 was relatively constant within basins of adjacent orders. 
We will now demonstrate that the constancy of this coeffi- 
cient is largely established by the constancy of the drainage 
(D). To do this we choose a,8 pairs representing the equi- 
librium solution within the a,8 admissable region. By 
feet ly Setting ay ee | vongescolying equation 3. | 


1 
together with the relationship 
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da = Lia ‘ (32) 


which is equation 1.6 with i = N and &. 


i = devs we obtain the 


values for ae and Ag set forth in Table 3-2. We are able 

to let ay and 2, equal unity here initially since we are 
interested in demonstrating the dependency of k on the drain- 
age density (D) and not in its absolute values. Consequent- 
ly, the values for Le and Ag set forth in Table 3-2 have no 


absolute significance. 


FABLE S= 2 


LENGTH AND AREAS OF AN EIGHTH ORDER BASIN 
FOR VARIOUS EQUILIBRIUM a,8 PAIRS 


a B g A 


8 8 
iO 64 422 
Brae 244 Eno 3S 
4.24 ve a, 18104 
ayes > 1838 82495 
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Now, by varying the drainage density we can calculate 
various sets of Le, Ag - values in the manner used to pro- 
duce those set forth in Table 3-2. We then perform a 


least-squares regression analysis on each set to examine the 


behavior of k and h in the power equation 
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The regression equations are 


2 lo Ag = F°L 16¢q Le rm LOO; 
Etleg Ag tog fo) ="F 5 (10g be)? + log b £ log 2g, 
and 
L 
pibae oF 
te = (2 AQ)” . 


The analyses give the results set forth in Table 3-3. 


TAREE 32:3 


RESULTS OF LEAST-SQUARES REGRESSION ANALYSES 
FOR VARIOUS DRAINAGE DENSITIES 


‘ By Changing L, By Changing ay 
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Thus, we see as did Hack [Ref. 27] that the coefficient (k) 
is a function of the drainage density (D) and that the 
exponent (h) 1s necessarily a function of the overland 

flow ratio (9). Actually, as pointed out by Hack, kK 1s alee 
affected by the valuesof p but to aiverys small cx@ecnt. New 
Since k is largely a function of drainage density we would 
expect k to remain essentially constant within a relatively 
homogeneous basin. In this respect, a number of authors 
have found that not only does k remain relatively constant 
within any given basin, but actually tends toward a common 
value of 1.4. In a study of streams in the Shenandoah 
Valley, Vaurggnaa, Hacks (Ref. 27] recorded the relatwonshie 


depicted by Figure 3-3. 
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Figure 3-3. Relation of Channel Length in Miles to Drainage 
Area in Square Miles (redrawn from Hack [Ref. 27]) 
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The relationship illustrated by Figure 3-3 has been verified 
by many authors in describing a number of widely scattered 
areas of the world. Thus, throughout the balance of this 


thesis we will refer to 1.4 as the equilibrium value of k. 


C. THE BASIN SHAPE PARAMETER 

Schumm confirmed the findings of most of the earlier 
authors in the field when he wrote "as the relief ratio 
increases the drainage basin becomes more elongate" [Ref. 
34: 23]. In describing the drainage network of the badlands 
at Perth Amboy, New Jersey, Schumm found that as drainage 
systems mature and approach an equilibrium state, the relief 
increases and causes the basin shapes to become more elon- 
gate. Since we can characterize a drainage baSin in its 
earliest stages of development as one in which no jumping 
has occurred, Schumm's findings imply that the process of 
maturing is closely allied with the growth of the jumping 
process. Thus, if we consider an Nth-order drainage basin 
in which jumping is prohibited, we would expect the basin 
shape parameter (h) to decrease toward 0.5 as the basin 
order decreases. The latter would be true since we have 
shown that h is unaffected by changes in the drainage 
density and since the jumping prohibition implies that all 
streams lie within the basins of the next higher order; the 
latter condition implying that both a and 8 remain constant 
within the included basins of orders less than N. Thus, we 
Surmise that as stream networks are in their initial stages 


of development, the typical first order stream basin is 
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Circular. As the network matures and larger streams are 
formed, the P. can no longer equal zero for all i<N and 

the basins begin to elongate in a downstream direction. 

This process continues until the network achieves a steady 
state and the basin shapes become relatively constant within 
the parent basin. This process can be illustrated by 
letting the elongation occur as a one-sided ellipse in the 
downstream Wie ction. The circular area and the one-sided 


elliptical area are illustrated in Figure 3-4. 


b}—— ly ——4J 


Figure 3-4. Typical Drainage Basin 


The area-length relationship of the typical basin can be 


expressed as 


ee Ty ba 
Ay = 4 Li (Li ae ( ea) 
which, with € = 0, is the equation of a circle. We can also 
write 
I Aal 
y eee 
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where 


— lo Er Fen (Sey = (3.9) 
Jo N Jo N 


Thus, relationship 3.8 becomes 
_ “e h 
Ly. = = Ay) ' (3.20) 


For a given length, h is a function of ¢€, the downstream 

elongation. As € increases from zero, h increases from 0.5 
towards unity but can never achieve unity since e<Ly- 
Mayer [Ref. 32] and Hack [Ref. 27] found that basin length 


can be more closely represented by 


i. =. | oo Lep 


which implies that the typical mature basin area is somewhat 
greater for a given length than is portrayed by the one- 
Sided ellipse. In any event, the area of the drainage 

basin does vary as the hth-root of the basin length, the 
value of h being a function of the overland flow ratio (p). 
Solving equation 3.11 for h and substituting therein 
equations 3.1 and 3.7 we can obtain the equation for h for 
the general case i. Then forming the ratio of adjacent 


basin shape parameters we further obtain 





- gt-et 
h. log. k& ee ee ee a- 8 
a ee ——————————————— (3.139 
ava — j-l i-l ,i-l1 ’ ; 
Jo 1 a-8 
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where k is the constant that converts the length of the meni 
stream given by equation 2.3 into the length of the basin 
given by equation 3.11, i.e., the main stream is k times 

as long as the basin. If we now consider a basin which has 
achieved equilibrium with its environment, we can see that, 
although a is unaffected, 8 decreases within included basins 
of successively lower order since streams that jump do not 
lie within the basins of the next higher order. We can 
describe the behavior of the bifurcation ratio-within the 
included basins by letting Bs represent the bifurcation 


ratio within the ith-order basin and observing that 
> oe 


where 8 is the constant bifurcation ratio for the parent 
basin and qi} the proportion of (1i-1)st-order streams that 


don't jump. Using equation 2.1 to obtain q._, we can write 


a (1-9) + 28 7 (p-pN7**4) - 2 (3.13) 
1 - Natl 

Now, if h is relatively constant for adjacent included 

basins as assumed in Chapter II, equation 3.12 should give 

solutions near unity as 1 varies when a remains constant 

and B. is defined by equation 3.13. Solution of the equation 

for N = 8 and the equilibrium values of a = 2.68, 8 = 4, 

and k = 1.4 disclose values of .99+ as N increases from 3 

to 8 indicating that the assumption is valid for included 


basins of order near N. Since equation 2.4 requires only 
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that h be constant for the Nth- and (N-lJst-order basins, we 
see that the requirement is fully met. This theoretical 
development is supported by a number of empirical investi- 
gations. The specific conclusion concerning the constancy 
of h has been substantiated by Hack who, in commenting on 
his equation L = bA’, stated 

the exponent n must always have approximately the 

same value in any homogeneous area for its value 

is determined by the laws of probability and not 

by geologic factors [Ref. 26:Bl7]. 
In examining streams in Michigan, Virginia, and Maryland, 
Hack [Ref. 27] further found that as the drainage basins 
became more homogeneous, h approached a value of 0.6 and 
k a value of 1.4. In their random walk model of a hypo- 
thetical drainage system, Leopold and Langbein [Ref. 28] 
found that a value of 0.64 held for h. This latter finding 
was Significant not only for the value of h but because it 
strongly supported Hack's finding that h was the product 
Of chance factors. In a Horton analysis of 47 small 
watersheds in Illinois, Iowa, Missouri, Nebraska, North 
Carolina, Ohio, and Wisconsin, Gray [Ref. 9] found that 
values of 1.4 and 0.568 held for k and h, respectively. 
Schenck [Ref. 18] used a Monte-Carlo random walk computer 
Simulation to show, among other things, that Gray's findings 
were substantiated by considerations of the most probable 
state. One of the more significant substantiations was 
contributed by Leopold, et al., [Ref. 1] by providing a plot 


of many of the world's major rivers showing an adherence to 
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a 0.6 value for h worldwide. Mayer [Ref. 32] observed that 
a constant h does not imply a constant length-to-width ratio 
and provided the values set forth in Table 3-4 as an 


A Liustration. 


TABLE 3-4 


LENGTH-TO-WIDTH RATIOS OF DRAINAGE BASINS 
(REDRAWN FROM MAYER [REF. 32]) 


h = 0.6 
Area Length Width Length/ 
(sq mi) (mi) (mi ) Width 

100,000 dle OORO) 100 iG rel 
107,00) 250 40 Ga2 
1 00u Ges 16 Sit 
100 16 Gre Pies: 
10 4 25D dyes 5 
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C. MEAN BASIN WIDTH 

Hack [Ref. 27] also observed that the width of a typical 
basin could be approximated by the'mean width (W) which, Li) 
effect, reduces the basin shape to a rectangle. Thus, we 


Can weice 


which substituted into equation 3.11 yields 
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Wy =o, 3 (3.14) 


Using the mean width would most likely yield unsatisfactory 
results if only one or only a few drainage basins were 
crossed. It would, of course, depend on where the basins 
were crossed. However, if a sufficiently large number of 
basins 1S considered and the points of crossing randomly 
determined, the expected crossing point would occur where 
the basin assumes its mean width. Thus, in this context, 
which is quite reasonable from a planning viewpoint if the 
standard deviation is not too great, the mean width is 


fully acceptable. 


Sil 


VI. A CROSSING ALGORITHM 


A. THE HORIZONITAL COMPONENT OF STREAM LENGTH 

We consider an Nth-order parent basin possessing mean 
characteristics in which Q. and n, again denote, respective- 
ly, the mean length of the ith-order streams and the number 
of ith-order streams within the mean basin area. We assuie 
that the streams are randomly distributed and randomly 
Oriented within the drainage basin and that 8, as given by 


Figure 4-1, is uniform on the interval (0, T/2). 


¥y 
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Figure 4-1. The Horizonital Component 
of Mean Stream Length 


We further envision a force of width W traversing the mean 
basin parallel to the minor axis and define the horizontal 
component of each mean stream length as that length perpen- 
dicular to the axis of traverse. The horizonital component 
is given by Q. cos @ in Figure 4-1. The density function 
of 8 can be written as 


727 ee Ds el yy. 
f(6) = ; (4.1) 


0 ; elsewhere 


ae 


The expected length of the horizonital component can be 


shown to be 


- 1/2 oo 
E[%. cos 6] = s |2, cos o| £(6) dé == &, (4.2) 


B. THE STREAM PROBABILITY DENSITY FUNCTION 

We denote the expected length of the horizonital compo- 
nent as te and consider the midpoint of each. Using the 
crossing force of width w we define the relationship por- 


trayed by Figure 4-2, where the coordinate system is 


envisioned to be superimposed over the mean basin area. 
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Figure 4-2. The Relationship Between Crossing Width and 
Midpoint of Horizonital Component. A midpoint to 
(fio wlohe wOoteZzenO Jandwto the right of wt; amplies 
that the respective stream is not crossed. The . 
midpoint location within the interval zero to wtj, 
inclusive, then, describes the length of stream 
protruding into the crossing zone given that the 
stream 1S crossed 
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Since streams are randomly distributed within the mean basin, 
the midpoints of he are uniform on the interval (O,w+h.). 

We denote the midpoint of an an as x. and define it as the 
position, measured from the origin, of the midpoint of the 
ith-order stream. We also define Y, as the extension into 
the crossing zone of the ith-order stream. We then let Xs 
and Y; represent the actual points and lengths, respectively, 


and define 
6, = set of indices i such that ds coetih. 
and O5 = setk®om iLndiseeshi eswch “that Le > Gn 


The value of the statistic va is then given by 
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The distribution of Key conditional on a stream being 


encountered, 1S given by 
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A < XA °< +R. 
W+h . 0 i. : 1 


g(x,)= (4.5) 


0 - elsewhere 


The expected length of an ith-order stream extending into 


the crossing zone is then given by 


3 - 
ue LW 
Ely. ] = ‘ h, (x, )g(x,) dx, = oF T, ; (4.6) 
ied, 
and 
uti, 7 
Ely, ] — er he (x, )g(x,)dx, os ck - (4.7) 
| 0 i 
1€0. 
The respective variances are given by 
ve ms 
QR. (2w-k.) 
5 = —_— ’ (4.8) 
Yi 3 (wth. ) 
1€9, 
and 
5 ia Oe 
O = re ee : (4.9) 
Yi 3 (wth, ) 
1€9, 


The foregoing could perhaps be made clearer by considering 


the density of ve given by 


aD) 


9, ly.) F (4.10) 


1€6 


and 
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Using these mixed densities it can be easily verified that 
Yi: defined over the two regions indicated, does, in fact, 
have the mean and variances indicated by equations 4.6 


ENLOuUGiiE 4s. 


G&. CONSTRUCTION OF THESMODEL 

We now define n,; as the number of ith-order streams 
extending into the crossing zone. The expected total 
length of ith-order streams extending into the crossing zone 


is then given by 


eves 

me a Etn.y,] = ~ = ° 
W +2. 
ae 





The behavior of stream lengths extending into the crossing 
zone as a function of midpoint position is illustrated by 


Figure 4-3. 
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Figure 4-3. Behavior of Stream Lengths in Crossing Zone 


Since, by assumption, the midpoints of streams are uniformly 
distributed within the mean basin area, we can consider that 


the actual density of midpoints per unit area is given by 


Now, if a midpoint is within the width w + Lay the corres- 
ponding stream extends into the crossing zone. Thus, the 
midpoint area traversed by the crossing force in crossing 


the mean basin along the minor axis iS given by 


Weealty + Sas 


and the number of ith-order streams extending into the cross- 


ing zone is, therefore, 


(4.12) 


es, 


The expected total stream length of ith-order streams 
crossed by the force in one traverse of the mean basin along 


the minor axis is, then, 


o° = AW Qe w (Aiken 


The foregoing development can perhaps be better motivated 
by considering the following argument. Assume the following 
as we have previously: 

1. » the fecal midpoints are uniformly distributed 
within the mean basin area. 

2. The crossing zone is randomly chosen in the mean 
basin area. 

3. In any portion of the crossing zone the crossing 
force encounters any stream having its midpoint within 
the width Ww + h and does not encounter any stream not 


having its midpoint within this width. 


Now divide the crossing zone into n equal portions indicated 
in Figure 4-4. If n is large enough so that most of the 
small areas are randomly related to any particular one, the 
chance of failing to encounter a particular stream during 
the traverse is the product of the chances that encounter 
fails during metionealéng pei small piece. Now the prob- 
ability that the stream is encountered in a small area is 
the probability that the stream is in the small area which 


is given by 
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Crossing 
zone 


=l 


Basin Area 


Figure 4-4. Division of Crossing Zone Into 
Arbitrarily Small Areas 


Hence, the chance that along the entire traverse the stream 


i NoLeeweountered is 


W R n 
ile We ate Le) ] 


Thus, as n gets large the probability of encountering the 


stream approaches 


s W(wtd, ) 


aye, 


Tis “fUrener COnvereseror Manor aAxUS Claverse to 


(Wh. ) 


Now, the latter two functions are both cumulative distri- 


butions for the exponential probability density functions 


; (wth, ) 
h 
Die (wth, ) a “N j 
and 
; (wth, ) 
1-h 
ae (wth, ) = *N 


Since the distance between encounters iS exponential with 


mean 





the underlying process must be Poisson with mean for a minor 


axis traverse of 


—h A 
Ay (Wh. ) ’ 
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and for a major axis traverse of 


h-1 ‘ 
Ay (wth. ) ; 


Hence, we can write the expected number of streams encoun- 


tered as 
ro Oa ~ 
a a Ay B (W+h, ) f 
and 
¥x alt-l 2Niem a 
ie _ Ay B (wth. ) t 


which, for the minor axis traverse, is precisely equation 
4.12. This development does not give us expected stream 
length in crossing zone but does give us the added feature 
of mean interfluvial distance, i.e., distance between 
streams, which, using the mean of the exponential distri- 
bution, can be expressed as 


AN 


qd, = oe eee ° (4.14) 
1 B > (ute) 


Now, from equation 1.1 we can deduce that 


However, this equation reflects Horton stream lengths which 
are unsuitable for our purposes. Converting them to 


Strahler stream lengths using equation 1.9, we get 
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which, uSing equations 2.3 and 3.11, further converts to 


— _ h oie ieee 
de = k Ay (a O ) raat See aes eee (45555 


Using the fact that le = 2 i 


A along with equations 1.5, 3.1] 


3.14, and 4.15 we can express equations 4.8, 4.9, 4.12, 4.13, 


and 4.14 in more tractable terms as 


3 


2k .h, i-N _i-N-l i-N-1 
3? _ la An’ = 6) Mm) eas AN (or OL )] 
-N yt Tie. 2 ae 
Yi 3 [w+ ~ AD (a = al 
1€0 (4 69 
ame 
ay i-N i-N-1l 
wf ar (art Nat Nt) - w) 
ee en ioe 2 da 
Yi 3[w + — Ay (a - a )] 
ies, 
x N-i -h 2k i-N i-N-l 
ns B [WAL + == (ao OL |e eee eee (4.0 153)) 
p= _ a” eee ee (4.19) 
and 
i-N Dio ahi eee LoN oleae 
d. = £ Ay [w + —— Ay (a = 10 Vy pale” (4.20) 
Note: For i= 1 in the foregoing equations set 


gr7N-Lt — 
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Equations 4.18 and 4.19 can also be made applicable to 
crossings of the mean basin along the major axis by sub- 
stituting basin length for basin width. The equations 


then become 


kk fl ae 2iow,Jda aL “ea Ib 


we @ Ay [w + —— Ay (a =—<¢ )] , (4.21) 
and 
os gN-i as 5 Nee (g27N a gi-N-1) (4,22) 


Two additional aspects of the model of interest are the 
coefficient of variation (CV, ) and an upper bound (U, ) on 
the total expected stream length (5). The upper bound can 
be established by noting that the expected extension into 
the zone of one stream cannot be greater than w or yy 
whichever is appropriate, and the equations for the upper 


Pound then written as 
imo We. LS , LESS , 
and 


ma nee ied. 


The coefficient of variation is the standard deviation 
divided by the mean and is given by the equations 
(2a ioear aes We 
Rs (2-2; ) 
CV; = a ; ied, i 


ote 
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and 





D. MODEL APPLICATION 

A number of methods are available to obtain values for 
the model parameters. If a, 8, h, and k are known for the 
particular area in question then the model equations can 
be solved merely by measuring the area (Ay) of the basin to 
be crossed from a topographic map or aerial photograph. The 
area may be determined by using a planimeter or estimated 


using the equation A, = L Wyr where W,, 1s taken as the 


N N N 


‘width at the intended point of crossing. However, the 
latter method is understandably cruder than using a plani- 
meter. In the event that either a or 8 is unknown, the 
other may be extracted from the appropriate curve in Figure 
2-2, that is, by solving equation 2.8 for the appropriate 
value of h. In the event that none of the parameters is 
known, results for the most probable state may be obtained 
by using h = 0.6, k = 1.4, and 8B = 4, and again solving 
equation 2.8 for the appropriate a-value. Using 8 = 4 

to represent the most probable state has been substantiated 
by a number of authors. Actually, B varies somewhat from 
basin to basin but clusters about 4 in mature basins. The 
sensitivity of the model with respect to the various 
parameters will be discussed in the next section. At this 


point two model applications will be examined. 
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1. Application 1. When planning for a combat force, 
Such as a division, to move across a drainage area the 
conditions throughout the entire crossing zone could 
easily be of interest. To obtain this information we 
solve the model equations letting w equal the actual width 
of the crossing force front. Solution of the model for a 
menmer axis traverse for w = 1 mile, h = 0.6, kK = 1.4, 
® = 2.68, 6 = 4, N = 8, and Ag = 1500 square miles is set 
forth in Table 4-1. The solution for a major axis traverse 
using the same parameter values is set forth in Table 4-2. 
In solving the equations of this model, it is important to 
note that R must be taken as min [2 Ww! for a major 
axis traverse and min ae Ly] for a minor axis traverse. 
This is necessary to preclude violation of the basin 
confines. 

For discussion purposes consider the entry in Table 4-1 
for a fourth-order stream. The model predicts that the 
eCighth-order basin will contain 256 fourth-order streams 
all of which lie entirely within the confines of the basin 
area. If the crossing force has a crossing front of one 
mile (w= 1), the model predicts that approximately six of 
the streams will intersect the crossing zone and that the 
total expected length within the zone perpendicular to 
the crossing axis will be 2.77 miles. Now, given that a 
fourth-order stream intersects the crossing zone, the model 
further predicts that the expected extension into the zone 


of each stream is 0.47 miles, or about one-half the width 
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of the crossing zone. The latter statistic implies that 
given sufficient maneuverability in zone it is conceivable 
that, if necessary, the fourth-order streams could be by- 
passed. A decreasing frontage (w<l) would serve to limit 
this possibility. The model then computes a standard 
deviation for the expected horizonital length of streams 

in zone*of 0.74 miles and an upper bound of 5.19 miles. 

The upper bound is associated with certainty, i.e., given 
that six fourth-order streams intersect the crossing zone 

we can be 100 percent certain that the total extension into 
the zone of all fourth-order streams does not exceed 5.19 
miles. For the particular probability density function used 
in this model (equations 4.10 and 4.11), the standard 
deviation for fourth-order streams is Le oteted with 
approximately 54 percent confidence, i.e., we can be 
approximately 54 percent certain that the total extension 

in Ee the. cles smig zone of all fourth-order streams lies 
between 2.03 and 3.51 miles. In fact, it can easily be 
shown that the maximum value attainable for the standard 
deviation is wo. /3 and that this maximum occurs when 

ies... An additional measure of the variation about the 
mean is provided by the coefficient of variation which is 
a dimensionless measure of the variation as a function of 


the mean. In general, the coefficient of Variation for 


this model behaves in limiting form pursuant to the bounds 
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However, the upper limit on the latter set of bounds 
envisions streams of infinite length which, of course, is 
rather impractical. Consequently, the upper limit would 
behave in a more restricted fashion as reflected by the 
results in Tables 4-1 and 4-2. In fact, the magnitude of 
the coefficients indicate that the mean is a rather good 
approximation of the distribution. Finally, the model 
estimates the distance between successive fourth-order 
streams as 3.13 miles. Now, the number of parent streams 
in the crossing zone (ng) has an additional and perhaps 
more enlightening interpretation. For the minor axis 
traverse of the eighth-order basin the number of parent 
streams is shown to be 0.57. This can be interpreted as 
the probability that the parent stream is crossed given a 
traverse at the point where the basin assumes its mean 
width. This is a particularly salient feature of the model 
Since it is usually the main stream that possesses the 
greatest obstacle value. 

2. Application 2. A combat force might be considered 
to move in one or more columms under circumstances that 
would require knowledge only of the streams intersecting 
the column(s), e.g., movement along a road or other confined 


route. This model accommodates these circumstances with 
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WwW = 0. Since length of streams in the crossing zone is of 


no concern under these circumstances only n., n 


., and d. 
i ak a 


are applicable. Solutions using the parameter values of 


application 1 for both axes are set forth in Table 4-3. 


TABLE 4-3 


MODEL SOLUTION FOR EIGHTH-ORDER BASIN FOR TRAVERSE 
ALONG BOTH AXES WITH CONFINED CROSSING FRONTAGE 


Stream No. No.Streams in Mean 
Order Streams Crossing Zone Inter- 
fia. ) in Basin +-— ——————_—___——_ fluvial 
(n.) Minor AXis Major Axis Distance, 
(ns) (ns ) miles 
(d.) 
1 16384 sig I I 63.54 re 
2 4096 6.18 26.69 5. On 
S 1024 yoga ig la'FS's 4.50 
4 256 27h dg es 6.72 
5 64 iF. 8.6 8.03 LOTUS 
6 16 INEZS 7.83 14.96 
7 4 Ooe oo u ya ggg 
8 1 0.56 IP lilt Soros 


fee oben ofNs finaly 

Since the model describes equilibrium behavior, that is, 
the behavior of a drainage network that has reached maturity 
or a steady state, it is instructive to examine the 


sensitivity of the model results with respect to the four 
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parameters upon which the model is based. A five-percent 
perturbation of the parameters will be used, however, it 
should be understood that, in nature, a change in one of 
these parameters 1S accompanied by changes among the others 
to maintain the system in equilibrium. Hence, as pointed 
out by Mayer [Ref. 32], a five-percent change one at a 

time imposes a condition of unnatural stress for the system. 
Table 4-4 presents the perturbed results for the fundamental 
model output. The equilibrium (unperturbed) solution is the 


last entry in the table. 


TABLE 4-4 


MODEL SENSITIVITY FOR THIRD-ORDER STREAMS IN 
EIGHTH-ORDER BASIN, MINOR AXIS TRAVERSE 





Variable Perturbed Results of Perturbation 
Perturbed Variable Third Order Streams 
Values 
* * 
we ihe = 
Bee 792 L305 22a 
8 
a2 1307 Bl 53 52 
255 1024 Ly eo See 
a 
2.81 ro24 16.08 3.36 
Ou ae, 1024 19.99 4.14 
h 
0.63 1024 n4.36 4.14 
v.33 ro7 4 16.66 3.93 
k 
1.47 1024 Lee 7 Awe 
Equilibrium 
Solution 
g 1024 16.87 Are 
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Obviously, the parameter to which the model is most sensi- 
tive is 8 as afive-percent change in 8 results in approxi- 
mately a 22-percent change in the model output. Consequently, 
the model output would be far more reliable if the B-value 
used was, in fact, the actual §8-value from the drainage 
basin being crossed. Another interesting feature of the 
sensitivity analysis is that although n; 1s sensitive toh, 
Oo: is not. This apparent anomaly is caused by the fact that 
the expected horizonital length (y.) of “the streams within 
the crossing zone varies in inverse proportion to the 
expected number (n5) of streams in the crossing zone and 
cancels out the effect. fThus, O° is not. a function of basin 
shape, a fact that is apparent from equation 4.19. The 
model would, of course, produce more reliable results for 

a specific basin if all parameter values were seeeetined 
from the basin being crossed. However, failing this, 
determining 8 and, perhaps, a for each basin would place the 


model results well within the realm of acceptability. 


F. CONCLUSIONS 

The model results should be viewed in full consideration 
of the assumptions underlying the model structure. A short 
review of those assumptions is instructive. In Chapter I 
we assumed that mean stream length, mean basin area, mean 
overland flow area, and numbers of streams are geometrically 
related to basin order, which, of course, is what makes the 
respective dimensionless ratios a, i, 9, and 8 constant. 


The validity of these geometric relationships is supported 


eZ 


by a great deal of empirical evidence and by considerations 
of the most probable state. However, the relationships 
were in large measure least squares approximations and 

some scatter did exist especially among the data from some 
of the less mature basins. We next assumed in Chapter II 
that streams bifurcate in direct proportion to the bifur- 
cation opportunity and defined bifurcation opportunity as 
the ratio of the total length of streams of each of the 
higher orders to the total length of streams of all of the 
higher orders in the parent basin. We were thus able to 
completely describe the jumping process in equation form. 
We then assumed constancy of the coefficient of drainage 
density (k) and basin shape parameter (h) in basins of 
adjacent order, assumptions that were later confirmed in 
Chapter III to a high degree of accuracy. These assumptions 
enabled us to define the interdependence among the dimension- 
less ratios a, 8, A, and h and to write a polynomial 
defining the admissable region for the a,f-pairs for each 
value of h. The assumption concerning bifurcation further 
enabled us to define the most probable a,f-pairs within 
each admissable region. In Chapter III we related the 
concepts of this paper to the works of a number of other 
authors and examined in detail the validity of many of the 
relationships employed. Finally, in Chapter IV we gave our 
model a statistical structure by assuming that streams are 
randomly located and randomly oriented within the parent 


basin and that stream encounters along either axis follow 


gs 


a Poisson probability law. The latter assumption simply 
implies that the number of encounters in a traverse is 
directly proportional to the distance travelled and that 
the distance between encounters follows an exponential prob- 
ability law. Using the mean conditions described by these 
probability laws we were then able to write equations in 
terms of our dimensionless ratios which yield such vital 
statistics as numbers of streams of each order crossed, 
total length of streams of each order protruding into the 
crosSing zone, distance between crossings, and to 
associate a confidence interval with each estimate. 

As indicated earlier the model is designed to be part 
of a larger planning model. The model obviates the need 
for detailed map and aerial reconnaissance in a specific 
scenario area involving a.large land mass and reduces a 
hitherto difficult problem to one that can be solved by 
a simple subroutine to the.main planning program. It is 
emphasized, however, that without parameter quantification 
in the specific area of interest the model merely describes 
the most probable state. In this connection Chorley wrote 
"Often the achievement of exact equilibrium in nature 
occurs only momentarily as variations about the mean take 
place, and in these instances the existence of the steady 
state can only be recognized statistically" [Ref. 25:137]. 
As indicated by the sensitivity analysis, especially with 
respect to 8, the variations about the equilibrium can be 


expected to be significant. Thus, it is important to 


74 


quantify as many of the parameters as possible in the area 
of interest. As a minimum, §8 should be so quantified, 
however, the model results can be made fully acceptable by 
quantifying both a and 8. Overprinting topographic maps 
with this information along with the order of the parent 
stream in the scenario area could be accomplished without 
great difficulty. A traverse across a scenario area would 
likely encounter a number of different a and 8 values. 
This would merely require that the model equations be 
solved for each a and 8 value for the specific traverse 
distance applicable to each. 

The model output would perhaps have its greatest 
utility when used in conjunction with a cost model and an 
effectiveness model in order to determine the best mix of 
stream crossing equipment for a unit of a given size oper- 
ating in a preselected scenario area. The model would 
facilitate comparison of inventory mixes and mixes of 
hypothetical equipment projected for future development. 
The model could even be used to develop doctrine and to 
select routes of advance within fixed scenarios. The latter 
aspects should render the model quite useful in the conduct 
of large scale war games, especially those involving compu- 
terization. For example, what would be the trade-off value 
in attaching a bridge company to Force A versus Force B 
when the two forces were advancing along different routes? 
A question of this nature is likely to occur within any 


war game, and for many scenarios could best be answered 


jus 


using the output of this model. The ability to answer 
questions of this nature without making weak assumptions 
and without the expense of actual stream counting should 
add a new dimension to war gaming, a dimension in which the 


entire war game would be more closely aligned with reality. 
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